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Comparison  of  Numerical  and  Experimental 
Time-Resolved  Near-Field  Hall  Thruster 

Plasma  Properties 

Ashley  E.  Gonzales,  Justin  W.  Koo,  and  William  A.  Hargus,  Jr. 


Abstract — Breathing  mode  oscillations  of  a  xenon  600  W 
Hall  effect  thruster  have  been  studied  using  temporally  resolved 
experimental  data  and  numerical  modeling.  Fluctuations  in  xenon 
neutral  near  infrared  (810-835  nm)  emission  in  the  near  field 
thruster  plume  have  been  measured  at  l-/is  resolution  using 
a  high  speed,  phase-matched  intensified  charge  coupled  device. 
Oscillations  in  electron  temperature,  3-9  eV,  have  been  inferred 
using  a  collisional-radiative  model  and  a  two-line  ratio  method. 
The  time-resolved  emission  and  electron  temperature  measure¬ 
ments  are  then  used  to  assess  the  accuracy  of  the  numerical  model 
HPHall.  Although  simulations  were  able  to  accurately  predict 
the  time  averaged  thruster  behavior,  the  model  greatly  under 
predicts  the  magnitude  of  the  oscillations.  General  phase  trends 
between  the  discharge  current  and  emission  as  well  as  electron 
temperature  are  consistent  with  observations,  suggesting  that  the 
model  is  capable  of  capturing  some  of  the  oscillatory  behavior 
despite  the  dampening  of  the  oscillations. 

Index  Terms — Breathing  mode,  hall  thruster,  HPHall,  plume 
emission. 

I.  Introduction 

HALL  thrusters  are  a  plasma  propulsion  technology 
widely  used  due  to  their  low  thrust,  high  specific  impulse 
operation.  With  increasing  demand  for  these  thrusters  to  per¬ 
form  a  wide  range  of  missions,  there  is  a  need  to  move  away 
from  the  costly  experiment-based  thruster  development  to  a 
more  efficient  numerical  modeling-based  approach.  Current 
efforts  have  focused  on  improvements  to  the  numerical  model 
HPHall.  The  HPHall  is  a  radial-axial  hybrid  particle-in-cell 
(PIC)  model,  which  is  based  on  a  fluid  treatment  of  electrons 
and  a  PIC  treatment  of  ions  and  neutrals  as  first  developed 
by  Fife  [1].  Although  HPHall  has  been  moderately  successful 
at  predicting  overall  thruster  performance  [1]— [7],  verification 
is  needed  to  determine  its  ability  to  accurately  model  the 
more  complex  features  of  thruster  operation.  Of  particular 
interest  is  the  low  frequency  (10-50  kHz)  fluctuation  in 
discharge  current,  often  referred  to  in  literature  as  the  breathing 
mode  [8],  [9].  The  fluctuations  are  believed  to  be  due  to 
a  periodic  cycling  of  the  neutral  and  plasma  density  in  the 
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exhaust  region  and  are  thought  to  have  a  strong  influence  on 
electron  transport. 

Although  simulations  have  been  able  to  qualitatively  repro¬ 
duce  the  plasma  oscillations  observed  experimentally  in  the 
breathing  mode  [1]— [4],  the  lack  of  time-resolved  data  has 
made  it  difficult  to  validate.  Recent  efforts  have  been  made 
to  develop  time-resolved  techniques  on  the  breathing  mode 
timescale  (30  jus). 

Recently,  time-resolved  probe  measurements  of  the  far 
plume  using  retarding  field  electrostatic  analyzers  as  well 
as  Langmuir  and  emissive  probes  have  provided  electron 
density,  electron  temperature,  plasma  potential,  and  ion  energy 
measurements  on  the  1— 10-// s  timescales  [10]— [13].  However, 
probe  measurements  are  limited  to  the  far-field  due  to  their 
interference  with  the  plume,  making  them  typically  outside 
HPHall’ s  computational  domain.  Optical  measurements  are  an 
attractive  alternative,  capable  of  nonintrusively  characterizing 
the  near  field,  with  resolution  on  the  ns  timescale.  These 
measurements  are  mainly  limited  by  signal-to-noise  ratio 
(SNR)  and  the  speed  of  the  optical  components.  Time-resolved 
emission  has  been  studied  using  2-D  CCD  images  of  the  plume 
as  well  as  interior  channel  measurements  using  a  thruster 
integrated  fiber  array  [10],  [14].  Additionally,  time-resolved 
laser-induced  fluorescence  (LIF)  techniques  have  been  devel¬ 
oped,  which  nonintrusively  measure  the  ion  energy  distribution 
function  [15],  [16]. 

In  this  paper,  time-resolved  emission  measurements  of  the 
near  field  of  a  600-W  Hall  thruster  are  presented.  Breathing 
mode  oscillations  in  xenon  neutral  near  infrared  (NIR)  emis¬ 
sion  are  measured  with  a  1-// s  resolution  using  a  fiber-coupled 
spectrometer  and  high  speed,  phase-locked  intensified  charge 
coupled  device  (ICCD).  Oscillations  in  electron  temperature 
are  inferred  using  a  collisional-radiative  model  (CRM)  [17]. 
The  resulting  characterization  of  the  breathing  mode  cycle 
is  then  used  to  assess  the  temporal  fidelity  of  the  numerical 
model  HPHall. 

II.  Experimental  Measurement 
A.  Test  Facility  and  Thruster 

The  measurements  reported  here  were  performed  in  vacuum 
Chamber  1  at  the  Air  Force  Research  Laboratory  at  Edwards 
Air  Force  Base.  This  vacuum  facility  consists  of  a  2.4-m 
diameter,  4.1-m  long  cylindrical,  nonmagnetic  stainless  steel 
vacuum  chamber  with  two  liquid  nitrogen  baffled  (76  K), 
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TABLE  I 

BHT-600  Hall  Thruster  at  Nominal  Xenon  Operating 
Conditions  and  Performance  [21] 


Parameter 

Value 

Anode  Flow  Rate 

2.45  mg/s 

Cathode  Flow  Rate 

0.197  mg/s 

Anode  Potential 

300  V 

Anode  Current 

2.05  A 

Magnetic  Current  (inner  coils) 

2.0  A 

Magnetic  Current  (outer  coils) 

2.0  A 

Thrust 

39  mN 

Specific  Impulse 

1530  s 

Propulsive  Efficiency 

49% 

1.2-m  flanged  gaseous  helium  two  stage  cryogenic  (15  K) 
vacuum  pumps  capable  of  pumping  speeds  of  48  500  L/s  on 
xenon.  A  cold  cathode  gauge  is  used  to  determine  a  chamber 
background  pressure  of  approximately  4  x  10-6  torr  (corrected 
for  xenon)  during  thruster  operation. 

The  Hall  thruster  used  in  this  paper  is  a  600  W 
BHT-600  with  a  3.2-mm  hollow  cathode  manufactured  by  the 
Busek  Company  (Natick,  MA).  This  thruster  has  been  studied 
previously  using  both  electrostatic  probes  and  various  opti¬ 
cal  diagnostics  [17]— [20].  The  BHT-600  has  an  acceleration 
channel  outer  radius  of  32  mm,  inner  radius  of  24  mm,  and 
depth  of  10  mm.  The  magnetic  field  is  produced  by  four  outer 
coils,  and  one  inner  magnetic  coil.  The  outer  and  inner  coil 
currents  are  independently  adjustable  for  optimization  of  the 
field  strength.  The  nominal  conditions  and  performance  of 
the  BHT-600,  as  provided  by  the  Busek  Company,  are  shown 
in  Table  I  [21].  Nominal  conditions  of  the  thruster  produce 
breathing  mode  oscillations  with  a  very  narrow  frequency 
(38  zb  0.3  kHz).  These  conditions  allow  for  the  phase-matched 
averaging  technique  described  in  subsequent  sections. 

B.  Emission  Measurements 

Plume  emission  measurements  of  the  near  field  were  col¬ 
lected  using  the  optical  system  shown  schematically  in  Fig.  1. 
Emission  was  collected  at  a  single  location  6-mm  down¬ 
stream  of  the  thruster  exit  plane  using  a  reflective  fiber  optic 
beam  coupler  with  a  6-mm  diameter  and  greater  than  97.5% 
reflectance  in  the  450  nm-2  jum  range.  The  emission  was 
then  transmitted  through  the  chamber  using  a  fiber  optic  feed 
through  and  a  200- ju m  fiber  with  a  numerical  aperture  of  0.22 
and  low  attenuation  in  the  NIR.  A  fiber  optic  adapter  was 
used  to  couple  with  a  Horiba  1250M  Series  II  spectrometer 
with  a  1.25-m  focal  length  (f/9),  300-G/mm  classically  ruled 
diffraction  grating,  blazed  at  500  nm,  and  a  dispersion  of 
2.6  nm/mm.  The  spectrometer  entrance  slit  was  fixed  at 
200  jum  to  increase  throughput  for  greater  signal  strength. 

The  emission  spectrum  is  detected  using  a  cooled  Andor 
iStar  ICCD  with  1024  x  1024  pixels  (13  jum  x  13  //m) 
capable  of  dark  currents  as  low  as  0.065  e-/pix/s.  The  camera 
intensifier  system  acts  as  a  fast  solid-state  shutter  capable  of 
gating  on  1-ns  time  scales.  Gating  is  controlled  by  transistor- 
transistor  logic  (TTL)  input  and  user  defined  gate  delay  and 


Passive 

Inductive 


Fig.  1,  Optical  system  used  to  measure  thruster  emission  using  a  1.25-m 
Horiba  spectrometer  with  attached  Andor  ICCD.  Emission  measurements  are 
triggered  using  a  custom  triggering  circuit  to  synchronize  measurements  to 
current  oscillations.  Full  breathing  mode  cycle  is  characterized  by  varying  the 
delay  between  the  trigger  and  ICCD  gating. 

gate  width.  The  individually  gated  shots  are  accumulated  over 
the  length  of  the  exposure  and  recorded  using  the  National 
Instruments  Lab  View  PC-based  data  acquisition  system. 

Emission  measurements  were  background  subtracted  and 
wavelength  calibrated  using  a  standard  xenon  spectral  lamp. 
A  relative  spectral  response  (RSR)  calibration  was  applied  to 
correct  for  any  wavelength  dependence  of  emission  detection 
due  to  the  transmission  losses  introduced  by  the  optical  com¬ 
ponents.  The  RSR  was  determined  using  a  National  Institute  of 
Standards  and  Technology  traceable  200-W  tungsten  filament 
standard  of  spectral  irradiance  and  comparing  the  measured 
gray  body  emission  with  the  ideal  spectral  response. 

C.  Timing  and  Triggering  Circuit 

Time  resolved  measurements  were  achieved  using  a  phase- 
locked  triggering  system,  which  synchronized  ICCD  gating 
with  1  jus  portions  of  the  breathing  mode  cycle.  Temporal 
resolution  was  achieved  by  varying  cycle  phase  of  the  gated 
portion.  Adequate  SNR  was  achieved  by  integrating  phase- 
locked  shots  over  a  0.5-s  exposure. 

A  custom  trigger  circuit  and  digital  delay  generator  (DDG) 
were  used  to  process  the  discharge  current  signal  into  an  ICCD 
trigger  compatible  TTL  output.  Thruster  current  oscillations 
were  measured  using  a  passive  inductive  probe,  which  acts  as 
a  bandpass  filter  (3-dB  points  of  120  Hz  and  20  MHz  with 
<1°  phase  shift  at  frequencies  between  1  and  100  kHz).  This 
low-level  voltage  signal  was  then  passed  through  a  low-pass 
filter  and  subsequently  differentiated,  so  that  the  zero  crossings 
of  the  differentiated  signal  corresponded  to  the  maximum 
and  minimum  in  the  thruster  anode  current  signal.  A  high¬ 
speed  comparator  was  then  used  as  a  zero  crossing  detector, 
which  produced  a  TTL  signal  synchronized  with  the  minimum 
of  the  discharge  current  signal.  The  TTL  was  then  passed 
through  a  DDG  to  output  a  clean  TTL  trigger  to  the  ICCD. 
Characterization  of  the  entire  cycle  was  then  achieved  by 
varying  the  DDG  delay  between  the  synchronized  TTL  input 
and  the  ICCD  trigger  output. 

To  ensure  gating  is  correctly  synchronized  to  the  desired 
portion  of  the  signal,  a  1-GS/s  oscilloscope  is  used  to  monitor 
the  thruster  discharge  current,  the  TTL  trigger,  and  the  ICCD 
gate  monitor  output.  A  sample  oscilloscope  trace  is  shown 
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Fig.  2.  Sample  oscilloscope  trace.  ICCD  settings:  1-/IS  gate  pulsewidth, 
13-//S  gate  delay.  Note:  only  80  jus  of  the  0.5-s  exposure  is  shown. 

in  Fig.  2.  The  discharge  current  is  shown  to  oscillate  at 
38  zb  0.3  kHz  over  a  mean  value  of  2.05  A.  The  sample 
trace  is  shown  with  ICCD  settings  of  1-//  s  pulsewidth  and 
13- jus  gate  delay.  Upon  receiving  the  input  trigger  and  waiting 
the  specified  delay  time,  the  ICCD  is  gated  for  the  duration 
of  the  pulsewidth.  Signal  is  accumulated  over  several  gating 
cycles  until  the  end  of  the  exposure.  Measurement  of  the  full 
breathing  mode  cycle  is  achieved  by  varying  the  gate  delay 
from  0  to  26  jus.  The  2  jus  delay  between  the  minimum  of  the 
current  signal  and  the  trigger  is  caused  by  the  delays  added  by 
the  active  circuit  components  and  was  found  to  be  consistent  in 
all  measurements.  Additional  delay  added  by  ICCD  electronics 
was  found  to  be  in  the  order  of  nanoseconds,  and  are  assumed 
negligible.  Time  zero  shown  in  subsequent  figures  will  be 
defined  at  the  fall  of  the  TTL  trigger,  corresponding  to  an  input 
gate  delay  of  zero.  A  ±0.2  jus  uncertainty  in  the  temporal 
resolution  will  be  used  based  on  the  half  width  half  max 
(HWHM)  of  the  current  FFT. 

D.  Electron  Temperature  Measurements 

Electron  temperature  measurements  were  achieved  using 
emission  lines  with  a  known  dependence  on  electron  temper¬ 
ature.  Emission  measurements  were  taken  in  the  810-840  nm 
NIR  region  to  include  the  823,  828,  and  834  nm  neutral  xenon 
lines.  A  collisional-radiative  model  (CRM)  based  on  the  model 
developed  by  Karabadzhak  et  al.  [22]  is  used  to  determine  the 
electron  temperature  dependence  of  the  hi3 /his  and  h34/his 
line  ratios.  The  Karabadzhak,  Chiu  and  Dressier  (KCD)  model 
uses  empirical  excitation  cross  sections  [23],  [24]  and  proba¬ 
bility  statistics  to  model  Xe  I  metastable  populations.  The  line 
intensity  per  unit  volume  can  be  approximated  by 

JL  =  ^l{N°Ne)BXCX  (1) 

where  Ni  are  the  species  number  densities,  and  the  rate  term 
Ba  can  be  determined  using 

=  ke0  +  - 2 — ^0  (2) 

where  a  =  N\/Ne,  and  kfj  is  the  excitation  rate  coefficient  at 
wavelength  X  for  the  i  —  j  species  collision.  The  metastable 
correction  term  CA  is  approximated  using  the  upper  2  pi  level 
degeneracy  and  branching  probabilities  [22]. 


Fig.  3.  Sample  emission  spectrum  of  a  BHT-600  at  nominal  operating 
conditions.  Measurement  taken  near  peak  of  discharge  current  oscillation  with 
1-//S  gated  shots  over  a  0.5-s  exposure  (19  000  cycles). 

Errors  with  this  temperature  method  are  due  to  two  main 
factors.  First,  experimental  uncertainties  in  the  emission  col¬ 
lection  will  propagate  through  to  uncertainties  in  electron  tem¬ 
perature.  Most  significant  of  these  is  the  uncertainties  in  line 
fit,  which  can  lead  to  line  ratio  uncertainties  as  high  as  5%  for 
low  SNR  line  profiles.  This  results  in  a  1  eV  uncertainty  for  the 
hi3 /his  ratio  method  and  a  2  eV  uncertainty  for  the  hi4 /his 
ratio  method.  Second,  there  are  some  uncertainties  inherent  to 
the  KCD  model  assumptions,  which  are  not  quantified.  These 
uncertainties  will  be  most  significant  for  the  hl3  line,  which 
relies  on  the  metastable  approximation. 

The  hl3 /his  ratio  method  relies  on  the  ratio  of  the  two 
strongest  of  the  lines,  reducing  the  temperature  uncertainty  due 
to  uncertainties  in  line  fit  associated  with  low  SNR.  However, 
this  method  relies  on  the  metastable  approximations  for  hl3 
line,  which  is  not  reflected  in  the  temperature  uncertainty. 
The  1%34/his  ratio  method  is  more  accurately  modeled,  but 
the  lower  intensity  of  the  h 34  leads  to  issues  with  SNR  and 
higher  line  fit  uncertainty.  Therefore  the  1%34/hls  ratio  method 
should  not  be  used  in  conditions  with  low  electron  temperature 
or  low  emission  signal. 

III.  Experimental  Result 
A.  Emission  Measurements 

Emission  measurements  were  taken  6-mm  downstream  of 
the  thruster  exit  with  a  temporal  resolution  of  1  ps  in  the 
NIR  (810-835  nm).  The  cooled  ICCD  achieved  a  high  SNR 
>24  dB  for  most  of  the  cases.  Emission  measurements  near 
the  minimum  of  the  discharge  current  were  difficult  due  to  the 
low  signal  (SNR  =1.9  dB).  A  sample  emission  spectrum  is 
shown  in  Fig.  3.  The  majority  of  the  measured  lines  in  this 
wavelength  region  are  strong  neutral  lines.  The  emission  from 
the  metastable  linked  823-nm  line  is  the  strongest  of  the  three 
lines.  The  metastable  state  allows  for  excited  neutrals  with 
long  lifetimes,  which  can  then  be  easily  reexcited  producing 
strong  emission  lines.  Additional  neutral  Xe  lines  are  observed 
at  820  and  827  nm,  but  their  low  signal  in  proximity  to 
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Fig.  4.  Normalized  Xe  I  line  intensities  in  comparison  with  discharge  current. 
Measurements  taken  6-mm  downstream  of  thruster  exit  plane  at  nominal 
operating  conditions.  Line  intensities  at  the  end  of  the  cycle  are  undetermined 
due  to  low  SNR. 

the  stronger  lines  makes  them  difficult  to  use  for  diagnostic 
purposes.  The  three  strongest  lines  (823,  828,  and  834  nm) 
are  the  focus  of  this  paper. 

The  peak  value  normalized  line  intensities  are  shown  in 
Fig.  4  with  a  representative  current  trace.  There  is  a  strong 
correlation  between  the  peak  of  the  discharge  current  and  the 
peak  of  the  Xe  neutral  emission,  as  observed  in  [10]  and  [14]. 
A2  /is  delay  is  shown  between  the  peak  of  the  normalized 
intensity  and  the  peak  in  the  discharge  current.  Based  on  the 
delay  and  travel  distance  to  the  anode,  an  8  zb  1-km/s  average 
axial  velocity  can  be  inferred.  Delays  measured  differ  from 
previous  measurements  by  Bouchoule  et  al.  due  to  differing 
thruster  geometry,  but  a  comparable  5  km/s  average  axial 
velocity  can  be  inferred  near  the  thruster  exit  using  a  similar 
method.  A  significant  drop  in  emission  signal  is  observed  after 
23  jus,  near  the  minimum  of  the  discharge  current.  The  low 
signal  in  this  portion  of  the  cycle  increased  the  uncertainty  in 
line  intensity  measurements,  especially  for  the  lowest  intensity 
line  at  834  nm. 

B.  Electron  Temperature 

Electron  temperatures  are  determined  using  /s23  /  ^828  and 
^834/^828  line  intensity  ratios  and  the  two-line  method  of 
the  KCD  model.  The  resulting  electron  temperatures  are 
shown  in  Fig.  5.  The  two  methods  are  in  general  agreement, 
with  electron  temperatures  fluctuating  between  3  and  10  eV, 
—4.5  ±1.5  jus  out  of  phase  with  discharge  current  oscillations. 
The  7823  / 7828  ratio  method  results  in  an  electron  temperature 
of  5. 6 ±1.0  eV  oscillating  at  ±1.7  eV  approximately  —3  jus  out 
of  phase  with  discharge  current.  The  7834/7828  ratio  method 
results  in  an  electron  temperature  of  6.3  ±  2.0  eV  oscillating 
at  ±2.3  eV,  approximately  —6  jus  out  of  phase  with  discharge 
current. 

Lower  temperature  measurements  of  the  hn/his  method 
may  be  due  higher  order  Xe  I  and  Xe  II  lines  being  integrated 
into  the  7823  signal,  causing  a  higher  ratio  and  lower  inferred 
temperature.  Future  measurements  will  look  into  increasing 
the  spectral  resolution  to  better  characterize  the  line  feature. 
Higher  scatter  in  the  h^/his  versus  the  hn/his  ratio 


Fig.  5.  Electron  temperature  fluctuations  6-mm  downstream  of  BHT-600 
thruster  exit  plane  at  nominal  operating  conditions  using  the  KCD  model  a 
two-line  emission  method. 

method  is  due  to  the  h^/his  method’s  higher  sensitivity  to 
uncertainty  in  line  ratio.  Low  signal  at  the  end  of  the  cycle 
(>22  jus )  led  to  higher  uncertainties  with  the  integrated  line 
intensity  for  all  three  lines,  and  therefore  resulted  in  a  higher 
uncertainty  for  both  temperature  methods. 

The  37%  ±  15%  oscillations  in  electron  temperature  are 
comparable  with  the  20%-30%  observed  in  far  field  probe 
measurements  [1 1]— [13].  However,  phase  (relative  to  7j)  vary 
across  studies.  Recently,  it  has  been  shown  that  the  magnitude 
and  phase  delay  of  the  electron  temperature  oscillations  may 
vary  immensely  depending  on  breathing  mode  conditions  [11]. 
Future  measurements  at  similar  breathing  mode  conditions  are 
needed  for  further  comparison  of  probe  and  optical  measure¬ 
ments. 

IV.  Numerical  Modeling 

HPHall  simulations  were  run  at  nominal  operating  condi¬ 
tions  with  a  time  step  of  0.25  jus  over  a  2-ms  time  period. 
A  sensitivity  analysis  was  performed  to  determine  the  model’s 
dependence  on  input  conditions,  particularly  sensitivity  to 
electron  mobility  profile.  HPHall  uses  an  effective  mobility 
profile  to  account  for  the  anomalous  electron  mobility  not 
described  by  classical  theory.  The  effective  mobility  can  be 
described  by 

=  y^_Lclassical  ±  y^Bohm  (3) 

Kb  1 

y^Bohm  —  “yjT  (4) 

where  Kb  is  Bohm  fitting  parameter,  typically  adjusted 
between  0  and  1.  The  Kb/  16  term  is  referred  to  as  the  inverse 
Hall  parameter.  HPHall  assumes  an  axially  varying  inverse 
Hall  parameter,  similar  to  that  measured  by  Meezan  et  al.  [25]. 

To  determine  the  model’s  sensitivity  electron  mobility  pro¬ 
files,  simulations  were  run  with  varying  inverse  Hall  parameter 
profiles,  as  shown  in  Fig.  6.  A  representative  test  case  was 
chosen  based  on  the  comparison  with  experimentally  observed 
performance  parameters  (7\  7sp,  /bm)-  The  results  are  then 
used  in  conjunction  with  the  KCD  CRM  to  predict  emission 
fluctuations.  Simulated  oscillations  in  emission  intensity  and 
electron  temperature  are  compared  with  the  experimental  data 
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TABLE  II 

Comparison  of  HPHall  Predicted  Performance 


Fig.  6.  Inverse  Hall  parameter  profiles  used  in  simulation  cases. 

to  assess  the  model’s  ability  to  characterize  the  breathing  mode 
behavior. 

V.  Numerical  Modeling  Result 
A.  HPHall  Simulation  Results 

A  summary  of  the  measured  and  simulated  thruster  per¬ 
formance  is  shown  in  Table  II.  Measured  thrust  and  /sp 
performance  metrics  are  based  on  published  Busek  values  [18]. 
Electron  temperature  predictions  from  HPHall  are  based  on 
an  electron  density  weighted  average  of  Te  evaluated  in  the 
experimental  emission  collection  volume.  Oscillation  magni¬ 
tudes  are  estimated  using  a  ratio  of  the  standard  deviation  to 
the  mean  value. 

Although  the  model  was  successful  in  predicting  the  overall 
thruster  performance  within  10%  for  all  four  cases,  all  cases 
had  difficulty  accurately  modeling  breathing  mode  behavior. 
Most  of  the  cases  were  able  to  predict  average  discharge 
current  within  10%,  but  all  were  significantly  under  the  32% 
oscillations  measured  in  discharge  current.  All  simulations 
predicted  higher  than  observed  electron  temperatures  with 
significantly  lower  oscillation  magnitudes,  2%-6%  versus 
37%  ±  15%  observed. 

Differences  in  the  four  cases  can  be  observed  in  the  breath¬ 
ing  mode  frequency,  /bm,  indicating  that  breathing  mode 
predictions  are  highly  sensitive  to  choice  of  inverse  Hall 
parameter  profile.  Changes  in  frequency  appear  to  be  non¬ 
linear,  suggesting  a  complex  relationship  between  the  electron 
mobility  and  breathing  mode  frequency.  Further  analysis  in 
this  paper  will  focus  on  Case  1 ,  based  on  the  compatibility  of 
measured  and  Case  1  predicted  values. 

Oscillations  in  plasma  properties  were  sampled  at  four  axial 
locations:  1)  channel  near  anode  (—5  mm);  2)  channel  exit 
(0  mm);  3)  near  plume  (6  mm);  and  4)  far  plume  (10  mm).  The 
cross  correlation  was  determined  by  evaluating  the  correlation 
coefficient  R  between  the  discharge  current  and  the  plasma 
property,  while  varying  the  delay  (r)  to  the  plasma  property 
as  defined  by 

„  ,  ,  TVd(t)-Id)(X(t  +  r)-X) 

Rld,x\T)  —  I - - - = —  (5) 

jTdd(t)  -  Id)2  Z(X(f  +  r)-  X)2 


Fig.  7.  HPHall  Case  1  cross  correlations  of  plasma  properties  with  discharge 
current. 


Cross-correlations  between  the  plasma  properties  with  the 
discharge  current  were  evaluated  at  each  axial  location  and 
are  shown  in  Fig.  7.  The  r,  which  corresponds  to  the  peak  in 
Rid,x  is  used  to  determine  the  delay  and  phase  shift  between 
the  discharge  current  and  plasma  property. 

A  summary  of  the  simulation  is  shown  in  Table  III.  The 
magnitudes  of  the  oscillations  are  estimated  using  ratio  of  the 
standard  deviation  to  the  mean  value  (x/x).  An  uncertainty 
of  zb 5  kHz  for  the  breathing  mode  frequency  is  based  on  the 
HWHM  of  the  spectrum  peak.  Delays  between  the  discharge 
current  and  the  plasma  property  were  determined  using  the 
peak  of  the  cross  correlation.  An  uncertainty  of  zbl  jus 
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TABLE  III 

Summary  of  HPHall  Case  1  Predicted 
Plasma  Property  Fluctuations 


Z 

(mm) 

X 

x/x 

(%) 

Jbm 

(kHz) 

1~delay 

(/iS) 

(deg) 

Id 

- 

2.1  A 

7 

34 

- 

- 

1.00 

-5 

9.4  1018  m’3 

4 

33 

-5.6 

-68 

0.75 

0 

2.2  1018  m’3 

5 

33 

-4.0 

-76 

0.26 

no 

6 

5.3  1017  m'3 

7 

28 

-6.4 

-65 

0.11 

10 

3.3  1017  m'3 

7 

33 

-4.9 

-59 

0.09 

-5 

l.ll  024  m  -'/s 

8 

33 

-1.1 

-14 

0.90 

0 

1.7  1023  m"3/s 

9 

34 

-0.4 

-5 

0.75 

TL'i 

6 

1.8  1022  m"3/s 

12 

34 

0.4 

5 

0.76 

10 

5.3  1021  m'3/s 

13 

33 

0.4 

5 

0.76 

-5 

1.1  1018  nr3 

7 

34 

0.4 

5 

0.91 

0 

6.3  1017  m’3 

8 

34 

0.4 

5 

0.93 

Tie. 

6 

3.4  1017  m’3 

8 

34 

1.9 

23 

0.90 

10 

2.4  1017  m’3 

9 

33 

1.9 

23 

0.89 

-5 

27.3  eV 

2 

34 

7.9 

96 

0.60 

0 

22.8  eV 

2 

34 

12.4 

153 

0.31 

1  e 

6 

16.4  eV 

4 

34 

-4.9 

-59 

0.48 

10 

12.4  eV 

5 

34 

-4.1 

-50 

0.58 

is  based  on  the  accuracy  of  the  cross  correlation  method. 
Phase  shift,  </>idiX,  is  determined  as  a  ratio  of  the  delay  to 
the  oscillation  period,  expressed  in  degrees.  Low  confidence 
values,  correlation  coefficient  ( R  <  0.30),  are  shown  in  blue. 

Oscillations  in  electron  density  and  ionization  rate  were 
found  to  be  relatively  in-phase,  \(j>\  <  23°,  with  discharge 
current.  Neutral  densities  were  found  to  be  out  of  phase 
with  discharge  current  oscillations.  This  is  consistent  with 
the  ionization  bursts  associated  with  breathing  mode  theory. 
Electron  temperatures  simulations  at  z  =  6  mm  were  found  to 
be  out  of  phase  with  discharge  current,  <p  =  —59°.  Although 
simulations  significantly  underpredict  the  magnitude  of  the 
electron  temperature  oscillations,  phase  behavior  is  compa¬ 
rable  with  phase  delays  observed  experimentally  between  the 
discharge  current  and  electron  temperature  (<p  =  —69°). 

B.  HPHall-KCD  Simulation  Results 

Mean  normalized  oscillations  of  emission  intensity  (Jx)  and 
discharge  current  (Id)  are  shown  in  Fig.  8  for  X  —  823  nm. 
Also  shown  are  oscillations  in  electron  and  neutral  populations 
(Ne,No),  the  rate  term  ( B ^),  and  the  metastable  correction 
term  (C^),  where  Jx  ~  NeN^B^C^ .  Oscillations  in  emission 
intensity  are  shown  to  be  primarily  due  to  strong  oscillations 
in  electron  density  with  additional  contributions  from  more 
subtle  oscillations  in  .  Oscillations  in  Bx  are  due  primarily 
to  oscillations  in  the  dominant  keo(Te)  term. 

Cross-correlations  analysis  shows  a  strong  correlation 
between  the  discharge  current  and  emission  intensity,  elec¬ 
tron  density,  and  the  rate  term.  A  weaker  negative  corre¬ 
lation  is  observed  between  the  discharge  current  and  the 
metastable  correction  term,  suggesting  the  metastable  pop¬ 
ulation  varies  inversely  with  discharge  current.  The  slight 
phase  delay  (</>  =  25°)  between  the  discharge  current  and 
emission  intensity  is  inconsistent  with  measurements,  where 
the  peak  in  intensity  occurs  before  the  peak  in  discharge 
current  (c j>  =  —28°). 


Fig.  8.  Predicted  oscillations  using  HPHall  coupled  with  the  KCD  CRM, 
with  emission  intensity,  Jx  ~  NeNoBACA . 


Fig.  9.  Comparison  of  measured  emission  intensity,  /^eas’  an^  HPHall-KCD 
predicted  emission  intensity,  T^dd- 

A  comparison  of  the  numerical  (HPHall-KCD)  and  experi¬ 
mental  emission  data,  normalized  for  Jx(Id  =  3  A),  is  shown 
in  Fig.  9.  Measurements  are  shown  for  one  breathing  mode 
cycle  due  to  the  phase-matched  experimental  procedure,  unlike 
numerical  data  which  is  shown  over  several  cycles.  Although 
simulations  greatly  under-predicted  the  magnitude  of  discharge 
current  oscillations,  the  HPHall-KCD  modeling  was  able  to 
predict  the  positive  correlation  observed  between  the  discharge 
current  and  emission  intensity.  The  shape  of  the  experimental 
curve  suggests  a  nonlinear  relationship  between  the  discharge 
current  and  emission  intensity,  emphasizing  the  importance  of 
cycle  phase  in  future  correlation  studies. 

VI.  Conclusion 

New  time-resolved  methods  have  been  developed  to  gain 
understanding  of  the  complex  dynamics  of  the  Hall  thruster 
breathing  mode.  Time  resolved  (1  jus)  emission  measurements 
were  taken  of  the  near  plume  using  a  high  speed,  gated  ICCD. 
The  7823/^828  and  /s34//828  line  intensity  ratios  were  used 
in  conjunction  with  the  KCD  CRM  to  determine  electron 
temperature.  Electron  temperatures  in  the  near  plume  were 
found  to  fluctuate  between  3  and  10  eV.  A  strong  correlation 
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between  the  discharge  current  and  electron  temperature  was 
observed  with  a  significant  phase  shift  (cp  =  —69°)  between 
the  cycles. 

Simulations  were  performed  using  HPHall  and  the  KCD 
model.  Results  were  compared  with  the  experimental  data  to 
validate  HPHall’ s  breathing  mode  predictions.  Although  the 
model  was  able  to  correctly  estimate  the  time  averaged  perfor¬ 
mance  metrics,  the  model  had  difficultly  capturing  the  more 
complex  breathing  mode  oscillations.  Differences  in  breath¬ 
ing  mode  frequency  between  the  simulation  cases  emphasize 
the  model’s  sensitivity  to  choice  of  inverse  Hall  parameter 
profile.  General  phase  trends  of  the  model  are  in  agreement 
with  the  experimental  data.  However,  the  simulations  greatly 
underpredict  the  magnitude  of  the  oscillations.  Additionally, 
the  simulations  predict  higher  than  observed  electron  temper¬ 
atures.  Similar  results  were  observed  by  Garrigues  et  al.  who 
compared  the  numerical  and  experimental  time-resolved  LIF 
measurements.  Although  the  model  is  able  to  capture  the  most 
probable  velocity  of  the  ion  velocity  distribution  function,  it 
underestimates  the  dispersion  of  the  distribution  by  as  much 
as  a  factor  of  four.  It  has  been  suggested  that  this  is  a  result 
of  the  imposed  time-independent  Bohm  mobility  profile. 

HPHall  attempts  to  capture  the  turbulent  mobility  using  an 
effective  mobility  with  the  addition  of  a  fixed  Bohm  mobil¬ 
ity  term.  The  Bohm  mobility  profiles  are  selected  through 
an  iterative  process  designed  to  capture  steady-state  thruster 
behavior  (including  the  discharge  potential  and  current),  but 
the  fixed  profiles  fail  to  capture  dynamic  thruster  behavior. 
In  reality,  the  turbulent  mobility  is  coupled  with  oscillating 
plasma  dynamics  and  therefore  some  interplay  is  expected  on 
the  10-100-kHz  timescale  between  the  plasma  density  and 
mobility  profile.  By  imposing  a  fixed  mobility,  it  is  effectively 
damping  the  oscillations  of  the  system,  which  in  the  case 
of  current  measurements,  may  also  be  producing  the  higher 
than  observed  electron  temperatures.  Future  modeling  work  is 
needed  to  develop  a  dynamic  mobility  profile,  which  efficiently 
captures  the  coupling  between  the  oscillatory  thruster  behavior 
and  the  effective  mobility  profile. 
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